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THE USE OF MOMENTS OF MOMENTUM TO ACCOUNT FOR CRYSTAL HABITS 

by 

Patrick G. Barber* 

Department of Natural Sciences, Longwood College 
Farmville, Virginia 23901, USA 

ABSTRACT 

A 3-step theory of crystal growth is proposed which involves first an 
association of molecules or ions in solution to form an impinging growth 
unit, then second the orientation of this unit prior to its impact on the 
surface of a crystal, and finally the attachment of this unit to the 
crystal face. From this theory the habit of a crystal is dependent upon 
the moments of momentum of the impinging growth unit. The results of 
sample calculations are presented for sodium chloride, succinic acid, 
sucrose, and chromium boride. The faces predicted by this proposed theory 
are compared with those predicted by other, energy-based calculations and 
with those experimentally observed. 

The proposed theory suggests alternative strategies for crystalli- 
zation and habit modification which may be of technological importance. 

Listings of the two computer programs that were used are provided in 
Appendix B. 


*0n temporary leave assigned to NASA Langley Research Center, Hampton, 
Virginia, 23665, USA, for 1984-1985. 



INTRODUCTION 


Crystal habits are economically important as well as being scientifi- 
cally useful for understanding the crystal 1 ization process. In commerce, 
many high volume solids such as sodium chloride and sucrose are sold in 
crystalline form. In electroplating, the success of a metallic deposit 
often depends upon the crystal form of that deposit. In the electronic 
high technology fields, the properties of semiconductors such as silicon 
and gallium arsenide are dependent upon crystal purity and habit. In 
medicine, the crystal habit and in vivo growth of crystals such as urea and 
weddell ite are important in diseases such as gout and kidney stones. Since 
the properties of materials often depend upon crystal structures, an 
understanding of the crystallization process itself and the various factors 
that influence this process are important. 

There have been a number of proposals made for understanding the 
crystallization process and the resulting crystal habit. An early, 
empirical method presented by J.D.H. Donnay and D. Harker (1) simply 
observes that crystal faces are usually the ones that have large 
interplanar spacings in the unit cells. The area of importance of a 
crystal face is proportional to the reciprocal of the interplanar spacing 
of each face as is summarized in relationship 1 below. 

A hkl a Vdfikl (1) 

Other more sophisticated methods of accounting for crystal habits have 
been developed based on Gibbs' principle that the crystal habit is the one 
with a minimum in surface energy (2,3). Typical of these is the molecular 
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theory of crystal growth which determines the importance of crystal faces 
from minima in free energies. Thermodynamics, molecular dynamics, Monte 
Carlo calculations, computer simulations, and graph theory (topology) with 
and without kinks, steps, and surface roughening have all been used to 
estimate crystal habits and growth (4-18). The Periodic Bond Chain (PBC) 
method of P. Hartman and W. G. Perdok (19,20) is typical of this type of 
calculation. Most of these methods assume that molecules are adsorbed from 
solution onto a crystal face and that the habit regulating event is the 
incorporation of the adsorbed species into the crystal. Such methods have 
been used to calculate equilibrium crystal habits in such compounds as 
sucrose, chromium boride, succinic acid, biphenyl, and weddellite (21-27). 
These methods focus attention on the bonding in a crystal or microcrystal- 
lite as the most important factor governing crystal growth. Much of what 
happens in crystal growth, however, occurs in solution during diffusion to 
the growing crystal. 


PROPOSED THEORY 


A 3-step, kinetically based theory of crystal growth is proposed to 
offer an alternative to the existing thermodynamically and statistically 
based ones. The first step is the association of molecules and ions in 
solution to form an impinging growth unit. The largest size of this unit 
is the crystallographic unit cell; but the actual impinging growth unit is 
expected to be smaller than this, since in most unit cells atoms and ions 
are shared among several adjacent unit cells. The formulation of this 
growth unit is a function of the degree of saturation, which in turn is 
dependent upon solvent, temperature, and ionic strength of the solution. 
Any proposed growth unit must, of course, be a chemically reasonable 
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association of ions or molecules that can geometrical ly build up the entire 
structure. Intermolecular forces can enhance the association process that 
originiates in the solution. This association step is assumed to be a 
rapid process which quickly attains equilibrium as indicated in figure 1 
and the following seven equations. 

Step #1: Association 


k l 

nA + mB <=> AnBm (random) (2) 

k 'l 

If the association quickly attains equilibrium, the previous 
equilibrium approximation can be made for the reaction in equation 2 above, 
and the following six equations can be obtained: 


k “lCA n B m (random)] = ki[A] n [B] m 
[A n B m (random)] = (kj/k-1) [A] n [B] m 

where 

kj = Aj exp (- A G^)/kT 
and 

k_i = A_i exp (- A Gti)/kT 


thus 

(k]yk_i) = (A 1 /A_ 1 ) exp [-( A - a G^i)/kT] 
From this it follows that 


(3) 

(4) 

(5) 

( 6 ) 


[A n B m (random)] = (A^A^) exp C(-E a )/kT] [A]"[BP (7) 

where 

E a = [ A - A G^] (8) 
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The concentration of the impinging growth unit is thus given by 
equation 7 and is a function of the solute concentration, the temperature, 
the activation energies, and other factors such as molecular geometry which 
are summarized in the pre-exponential factors, and A.j. The unit 
formed, A n B m , is randomly oriented in the solution unless the geometry or 
conditions of crystallization impose a preferred orientation. Such an 
imposition may occur in gel growth and in MOCVD systems which have 
streaming reactant solutions in which solute concentrations are 
inhomogeneous. 

There is some experimental evidence that these impinging units exist 
in solution. The existence of large molecular associations is accepted for 
colloids, and there is evidence that aggregations occur on smaller scales 
as reported by J. J. Burton (28,29), M. Hoare and P. Pal (30-32), and D. 
Bonchev, 0. Mekenyan, and H. Fritsche (33). 

Two simplifying assumptions are made at this point in the proposed 
theory. It is assumed that growth units do not change during or after 
attachment to the growing crystal face, and as such they have the same 
molecular arrangement in solution immediately prior to impact that they 
have within the crystal after incorporation. It is further assumed that 
only one type of impinging growth unit occurs in a given solution. 

Step #2: Orientation 

This randomly oriented unit must migrate from the bulk of the solution 
to the face of the growing crystal; and while migrating, this unit must 
acquire the correct orientation so that upon impact with the crystal it can 
attach. This second, orienting step is assumed to be the rate determining 
one. 
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k 2 

A n B m (random) < = > A n B m (oriented) (9) 

k_2 

The following equations are a consequence of the above reaction step: 
k 2 [A n B m (random)] = k_ 2 [A n B m (oriented)] hk i (10) 

C A n B m (oriented) ]hkl = (k 2 /k_ 2 ) [A n B m (random)] (11) 

[A n B m (oriented )] hk i = (k 2 /k_ 2 ) ( A i/A_i) exp (-E a /kT) EAHBJ" (12) 

The rate of growth-unit disorientation depends upon the moments of 
momentum of the migrating, impinging growth unit. It is assumed that 

k -2 = C / J hkl ( 13 ) 

where c is a constant of proportionality and J^i is the moment of 
momentum of the impinging growth unit about the axis through the center of 
mass perpendicular to the hkl-plane. The larger this moment of momentum 
the less easily the impinging growth unit will be moved away from the 
hkl-orientation needed for the proper impact, and thus the faster will be 
the growth of the corresponding face. Each hkl-direction and each 
molecular association has a different rate constant for disorientation, 
k_ 2 . In equation 13, however, it is assumed that this difference in rate 
constants is due only to the difference in moments of momentum about each 
hkl-axis. Thus, the constant of proportionality, c, in this equation is 
the same for each orientation. As a consequence, equation 12 can be 
rewritten. 

[A n B m (oriented)] hk i = (k 2 /c) (A^A.i) exp ( -E a /kT) (J hki ) [A]"^] 1 " (14) 
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From equation 14, the concentration of the correctly oriented, 
impinging growth unit is dependent upon its angular momentum. When this 
moment of momentum is small, the concentration of the impinging growth unit 
having the correct orientation for subsequent attachment and growth also 
will be small. Crystal faces that grow slowly are the ones that are 
observed ultimately to be the main facets of the crystal habit, since fast 
growing faces grow themselves out of existence. Thus, axes having the 
smallest moments of momentum in the impinging growth unit will be the same 
ones observed as facets in the resulting crystal. 

Step #3: Attachment 

The third step of this theory is the actual attachment of the 
correctly oriented impinging growth unit onto the growing face. 

*3 

A n B m (oriented) <— > A n B m (attached) (15) 

k-3 

The energies of the intermolecular bonds formed in the crystal are 
larger than the association or orientation energies of the growth unit in 
solution; and the rates of this third step have been assumed to be faster 
than the preceeding two steps, since the second step was considered to be 
the rate determining one. During crystal growth, a crystal can grow only 
as fast as correctly oriented material arrives at the growing surface, and 
all material is incorporated into a crystal upon reaching a growing 
surface. During dissolution of a crystal, however, the situation is 
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different. The energy of attachment in this third step may be the rate 
determining one during dissolutions, since the bonds holding the units into 
the crystal must first be broken before material can be carried away from 
the surface. 

As a consequence of reaction 15, the following two kinetic equations 
can be obtained. 

[A n B m (attached)]^! = (k 3 /k_ 3 ) [A n B m (oriented)] hk i (16) 

[A n B m (attached)] hkl = (k 3 /k. 3 )(k 2 /c)(A 1 /A. 1 ) exp(-E a /kT) (J hkl )[A] n [B] m (17) 

The ratio of the rate constants, k 3 /k_ 3 , for the attachment step is 
dependent upon the strength and energy of that attachment, which is a 
function of the number and strength of the bonds holding the unit to the 
growing face. This ratio is a function of the particular hkl-face being 
considered. Although ratios of rates for the forward and reverse reactions 
of different faces are not necessarily equal, it greatly simplifies the 
subsequent calculations if one assumes that these k 3 /k _3 ratios for the 
different hkl-faces are approximately equal. Because of this approximation 
this theory will be used to order crystal faces only among those faces 
which the Periodic Bond Chain (PBC), Burton, Cabrera, and Frank (BCF), or 
other energy-based theories have indicated are faces in which there are 
strong intermolecular forces of attraction to hold the added units to the 
crystal and are thus faces for which one may assume that the orientation 
step is rate determining. If all possible impinging unit orientations were 
considered, rather than just limiting the calculations to those 
orientations with favorable energies for attachment, one could easily 
obtain misleading results. Some high order planes may be found which have 
favorable moments of momenta indicating that the concentration of these 
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orientations are small and that the corresponding faces should be dominant 
in the crystal. However, it does little good to find planes of favorable 
momenta if upon impacting on the crystal surface the energies do not favor 
attachment. In such cases the principal assumption of this proposed theory 
has been violated, and the orientation step (step #2) is no longer the rate 
determining one. 

The rate of growth of a face, R^kl* 1S proportional to the 
concentration of the attached, correctly oriented impinging growth units. 

R hkl = k 4 E A n B m (attached)] hkl (18) 

R hkl - («4) ( k 3/«-3) ( k 2/lO (Ai/A_i) exp (-E a /kT) (J hkl ) [A] n [B] m (19) 
R hki = K exp (-E a /kT) (J hk i) [A] n [B] m (20) 

The above three equations have a constant of proportionality, k^ 
which is needed primarily to convert from units suitable for concentrations 
to those suitable for rate of growth. The differences in rates of 
attachment among the various hkl-faces have been considered above in the 
discussion of step #3. Consequently, k 4 is not dependent upon the 
hkl-di rection being considered. 

Since slowly growing crystal faces are the ones that ultimately have 
the greatest areas, the area of a particular face, A^i, is proportional to 

R hk l -1 » and 

A hkl/ A hkl' = R hkl / R hkl (21) 

where the primed factors are those for the most prominent hkl-face. 

Thus, by equation 21, the ratio of areas of crystal faces should be 
proportional to the inverse ratio of the rate of growth and thus inversely 
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proportional to the moments of momentum about the axes in the impinging 
growth units. A ratio parameter, r hkl , is defined in equation 22. 

r hkl = A hk i/A hk i' = Jf, kl /J hkl (22) 

SAMPLE CALCULATIONS 


By equation 22 the relative surface areas of crystal faces could be 
calculated from the moments of momentum of the impinging growth units. The 
calculation involves a computation of the moments of momentum about axes 
for several expected impinging growth units. These growth units should be 
reasonable molecular associations derived from the crystallographic unit 
cell. The same set of possible crystal faces considered by other theories 
will be considered in this paper. The moments of momentum will be used to 
rank the planes in order of expected importance for each model. A summary 
table is presented for each of the four compounds used as examples in this 
paper. The tables give the J-values, r-values, the order of the most 
prominent faces predicted by other theories of crystal growth as well as 
the experimentally observed order. 

The calculations were performed using two programs written in BASIC 
and run on an IBM-PC computer. Listings of the programs is given in 
Appendix B. The first program calculates the direction cosines of the 
vectors through the center of mass about which the moments of inertia are 
computed for each hkl-plane of the assumed growth unit, and the second 
program calculates the moments of momentum using equation 23. A derivation 
of equation 23 is given in Appendix A. 

J hkl - (MU) 2?=l(M I )[(x 1 ) 2 + (2/5)(ri) 2 ] (23) 
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Since literature summaries of observed crystal habits list only the 
most prominent crystal faces in a qualitative manner, quantitative analysis 
was not possible. The tables summarize the experimentally observed order 
of facet prominence as reported in the literature. Plus marks are used in 
the tables to indicate the extent of predominance if such information is 
available from the literature. 

The calculated values of the moments of momentum were obtained from 
equation 23 with the value of the multiplier factors, MU, selected to scale 
each calculation to the same number of atoms or ions as found in the unit 
cell. The values of were calculated by dividing each moment of 

momentum into the smallest value calculated for that trial impinging unit, 
as indicated by equation 22. 

Sodium Chloride, NaCl: 

Table 1 summarizes the results for the calculations of the moments of 
momentum for several impinging growth units for the sodium chloride 
structure. The impinging growth units are illustrated and labeled in 
figure 2. The crystal structure is the rock salt one with a = 0.563 nm 
(34). 

From this table and accompanying figure, the impinging unit that 
likely associates in solution prior to impact and growth consists of the 
simple association of Na + and Cl" ions into (NaCl)j units. As the degree 
of association increases (100) and (111) are more likely faces to appear. 
The associated growth unit consisting of 2 Na + and 2 Cl" ions arranged in a 
plane as (NaCl)2 is the structure used for model B in table 1 and 
illustrated in figure 2. Model D consists of exactly half the unit cell 
with an even greater degree of association and favors the (111) faces. 
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The moments of momentum for a few higher order faces were also 
calculated in Table 1. These terms all have large moments of momentum, 
which is in agreement with the observation that higher order faces are not 
common in crystals. 

A. Julig and B. Deprick have observed that the (110) faces appear and 
become more predominant as the degree of supersaturation increases in 
aqueous sodium chloride solutions containing glycine (35). As the degree 
of supersaturation increases one would anticipate a greater degree of 
association, and this theory predicts that (110) faces should become more 
likely, as the degree of association increases from model A to model B in 
table 1. Julig and Deprick also observed that (110) faces become less 
likely at any degree of supersaturation as the pH of the solution deviated 
from the isoelectric point of glycine. Buffers and electrostatic charges 
in solution can alter the degree of association and thus the resulting 
habit. 

Sodium chloride has been subjected to extensive studies to ascertain 
the effects of various additives on the resulting crystal habits (36). 

These studies conclude that selective adsorption of the additive onto 
crystal faces is the modus operandi . This could still be the mechanism 
under the theory proposed in this paper; but the proposed theory would 
allow for two further mechanisms for crystal habit modification. Altered 
habits could result from the attachment of the additive to the associated 
impinging growth unit while it is in solution. Such an attachment would 
alter the moments of inertia of this unit; or, the additive substance could 
alter the solution properties enough to cause a change in the size and 
shape of the impinging growth unit thus also altering its moments of 
momentum and thereby the resulting crystal habit. 
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The work of J. D. Birchall and R. J. Davey (36) on the effects of 
polysaccharides on the crystallization of sodium chloride suggests that the 
polymer adsorbs on the crystal surfaces and on the heterogeneous nucleation 
sites in solution. The theory presented in this paper suggests that the 
role of polymers in altering the rheological properties of solution should 
be considered also, especially since they conclude that "the polysac- 
charides are most active in influencing nucleation and growth when in chain 
extended configurations." This is just the configuration in which polymer 
solutions have the most altered rheological properties. 

Succinic Acid, C 4 H 5 O 4 : 

J. S. Broadley, D. W. J. Cruickshank, J. D. Morrison, J. M. Robertson 
(FRS), and H. M. M. Shearer reported the refined crystal structure of 
B-succinic acid in 1959 (37). Using their data, the moments of momentum 
and the r hkl“ values for &-succinic acid were calculated by the programs 
described in this report. The results are summarized in Table 2 , and the 
various impinging growth units considered are illustrated in Figure 3. 
Impinging units B, C, and E, consisting of molecules oriented along the 
c-axis, agree with the experimental observation that ( 001 ) is the most 
prominent face (24). These are reasonable molecular associations, since 
the c-axis is also the direction for hydrogen bonding in succinic acid. 
Model D has (010) predicted as its most favored face. 

The PBC analysis done according to the calculation method of Hartman 
and Perdock by R. J. Davey, J. W. Mullin, and M. J. L. Whiting (24) 
indicated that ( 110 ) should be the most pronounced face, but as they noted 
this is not in agreement with the observation that ( 001 ) is observed to be 
the most prominent face. The proposed theory of this paper using model C 
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predicts that (001) should be the face with the largest area and that the 
next face should be the (1,1,-1). The third most favored face in this 
model should be (110). This ordering of faces is in better agreement with 
the experimentally observed order of (001)>>(010)>(11-1)>(110) than is the 
ordering predicted by the PBC method. 

Impinging unit B emphasizes the ionic intermolecular forces at the 
ends of the molecule parallel to the c-axis. A review of the gamma values 
for this model in Table 2 indicate that this model may overemphasize the 
(001) face in the anticipated crystal habit. Nevertheless, Davey, Mullin, 
and Whiting (24) observed that the crystals of succinic acid become 
needle-like when grown out of isopropanol. This solvent has a smaller 
dipole moment than water, and perhaps the theory of this paper would 
suggest that isopropanol encourages the formation of another, competing 
impinging growth unit such as B, which is expected to favor more greatly 
the formation of the (001) face. 

Sucrose, Ci2^22^11 : 

Using the crystal structure of sucrose reported by George M. Brown and 
Henri A. Levy (38), the values for Jhkl and were computed and 

summarized in Table 3 and figure 4. 

0. Aquilano, M. Franchini -Angela, M. Rubbo, G. Mantovani, and G. 
Vaccari (21) observed the disagreement between the rates of growth 
predicted by a PBC analysis which they did using the methods of Hartman and 
Perdok and the experimental data for sucrose. From Table 3, impinging 
growth unit B agrees with the reported observation that the prominent face 
of sucrose is the (001) face. It is interesting to note that both the PBC 
calculation of Aquilano et al. and the BCF calculation of Saska and Myerson 
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assign a low value to the (001) face, yet this is an experimentally 
observed one. Saska and Myerson conclude that "The agreement of the 
calculated growth form with observed habits is fairly good and the 
reduction of the growth velocity in the c-axis direction [which is likely 
due to the preferential water adsorption on the (001) face] would result in 
an excellent agreement with the observed forms." [Page 554 of (22)]. The 
moment of momentum calculations predict the most prominent face and some of 
the other experimentally observed ones for sucrose without recourse to 
adsorption of water onto one of the faces. 

Chromium Boride, CrB: 

Using the data cited by R. Hamar and S. Hamar-Thibault (23), the 
moments of momentum for the faces for which they calculated the energies 
with a PBC calculation are summarized for several impinging growth units in 
Table 4. Various impinging growth units considered in these calculations 
are illustrated in Figure 5. 

From this table, impinging growth unit B agrees best with the observed 
crystal structure predicting not only that the faces (010) and (130) should 
be most pronounced but that they should both be about equal in surface 
area. This seems to be the observed order of preferred faces as reported 
by Hamar and Harmar-Thibault (23). Unit B in figure 5 has the two chromium 
and the two boron atoms arranged in a pseudotetrahedron. 

DISCUSSION AND CONCLUSIONS 


A 3-step theory of crystal growth has been presented in this report. 
The steps of association, orientation, and attachment occur in that order 
during crystallization with the orientation step being the rate determining 
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one. From this, a working hypothesis results by which the relative areas 
of facets of crystals can be predicted from the calculated moments of 
momentum of the associated, impinging growth unit using the I^kl of 
equation 22. 

The qualitative results of the moment of momentum calculations for 
several impinging growth units in the crystal structures of sodium 
chloride, succinic acid, sucrose, and chromium boride are compared with the 
observed order of importance of crystal faces. From this comparison, it is 
possible to select a chemically reasonable impinging growth unit that 
agrees with the observations. This theory allows for events occurring in 
solution to contribute to regulating crystal growth. Crystal habits may 
not be only a consequence of the thermodynamics of the solid. 

This proposed theory also provides alternative explanations for the 
mode of action of modifiers of crystal habits. Additive substances can 
adsorb onto the face of growing crystals as advocated by existing theories 
of crystal growth, but two additional modes are possible in the theory of 
this paper. The additives can adsorb onto the associated, impinging growth 
unit while it is still in solution and in this way modify the moments of 
momenta of this unit and thereby the resulting crystal habit. The 
additives can also modify the properties of the solution such as its 
viscosic or electrolytic properties to induce the formation of an altered 
impinging unit which will have different moments of momenta and an altered 
crystal habit. 

The results of calculations on four crystals have been summarized in 
tables and figures. The results of this theory are as good as the 
predictions made by other methods. The method easily allows for the 
possibility of quantitative comparisons and more detailed refinements of 
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the structures of impinging growth-units when relative crystal face areas 
are available from experiments. 

The modified model of crystal growth encouraged by this theory also 
suggests technological improvements in methods of crystal growth. Crystals 
of preferred orientations may be obtained if one is able to encourage the 
formation of appropriate associated, impinging growth-units in solution 
prior to crystallization. Such considerations may affect the growth of 
simple and compound semiconductors allowing desired crystal orientations to 
be more easily obtained. For example, large-area crystalline silicon from 
the catalytic decomposition of silicon compounds may be more easily 
obtained by selecting the decomposing compound or the gas flow conditions 
to encourage the formation of appropriate impinging growth-units. Similar 
considerations may assist in the growth of single crystals of GaAs in MOCVD 
systems and in epitaxial semiconductor growth studies. 

The theory as presented assumes the existence of a unique impinging 
growth unit. In solutions there may be a statistical distribution of 
several impinging growth units. Also, the energies of attachment may alter 
the molecular arrangements within the growth units upon attachment. Thus, 
the molecules in the impinging growth unit need not have the same molecular 
arrangement as the final solid. Consideration of these two additional 
factors may improve quantitative agreement between predicted and observed 
crystal habits. 
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TABLE OF SYMBOLS 


[A] and [B] 
A i 

A hkl 


d hkl 

E a 

AG^ and AG^i 
J hkl 

k 

k i 

k.i 

K 

Mi 

MU 

n 

R hkl 

T 

x i 

r hkl 


molar concentration of species A and B 

pre-exponential factor in the rate constant expression 

surface area of the hkl-face in a crystal 

a constant of proportionality between the rate constant for 
disorientation and the moment of momentum 

distance between hkl-planes in the crystal unit cell 

activation energy for formation of impinging unit 

free energy of activation for forward and reverse chemical 
reactions respectively 

moment of momentum for molecular model rotated about an 
axis through the center of mass and perpendicular to the 
hkl -face 

Boltzman's constant 

kinetic rate constants for forward chemical reactions 
kinetic rate constants for reverse chemical reactions 
a constant equal to (k4)(k<j/k-3)( k 2/ c )( A i/ A -i) 
mass of atom i 

normalization factor for each model 
covalent or ionic radius of atom i 
rate of growth of the hkl-face 
absolute temperature 

perpendicular distance from the hkl axis of rotation to 
atom i 

ratio of the area of an hkl-face to the area of the largest 
or most prominent face. 
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APPENDIX 1: The Calculation of the Moment of Momentum 


The moment of momentum for a collection of n particles is defined in 
physics texts (40) to be 


n _ __ 

J = £ x-j x m-jV-j (Al) 

i =1 

where x-j is the perpendi cular vector from the i-th particle to the axis 
of rotation. 

Since 7-j = to x Xj , (A2) 

it follows that 

J = Z-j[mj "x-j x (oj x x-j)] (A3) 

which contains a triple cross product made up of moments and products of 
inertia, which are terms having the forms given below: 

^i,x = E m i (yi 2 + z i 2 ) (A4) 

i 

and 

P X y = Ei m-j x-j y-j 

In the case of an associated impinging unit, it is assumed to be a 
rigid collection of spherical atoms, each one of which has a uniform 
density. The contribution that each spherical atom makes to the total 
moment is small compared to that contributed by the rotation of these same 
atoms about the axis of rotation. This is because the distance from the 
axis of rotation, x-j, is large compared to the radius of each atom, r j ; and 
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these distances enter the calculation squared. The contribution of each 
individual atom to the moment is given by 

(A5) 2/5 m^- 2 

where m^ is the atomic mass of the i-th atom (40). 

The contribution of each atom to the moment is the sum of its atomic 
component plus that arising by virtue of the rotation of this atom about 
the axis of rotation normal to the hkl -plane. The total moment for each 
atom is thus 

(A6) l,i - m i x i 2 + 2/5m-j ^i^* 

where r^ is the atomic radius and x-j is the perpendicular distance from 
the i-th atom to the axis of rotation. 

The total moment for the hkl -plane is the sum over all the atoms in 
the model being considered, or 

(A7 ) Jftkl - Z-jm-j [x-j 2 + 2/5r-j 2 ] 

Since each model may have a different number of atoms, it is 
convenient to introduce a scale factor, MU, which scales the total moments 
to the number of atoms found in the crystallographic unit cell. This does 
not change the relative order or crystal faces within a model but does 
facilitate comparisons among models. As an example, if a model has four 
molecules while the crystallographic unit cell has only two, then the scale 
factor MU is 1/2. Thus, 
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(A8) 


Jhkl ■ (MU) Ei-llMjlUx,-) 2 + (2/5)(r,)2] 


Equation A8 is the same as equation 23 used in this report. 
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APPENDIX 2: Computer Program Listings Used for the Calculation 

of the Moments of Momenta 

This appendix contains the listings of the computer programs used to 
calculate the moments of momenta described in this report. The first 
program computes the direction cosines to the normals of the selected 
crystallographic planes, and it computes the areas of those planes. The 
flow diagram for this program is given in figure 6. 

The second program listed in this appendix uses the direction cosines 
computed for each crystal lographic plane and the atomic positions of the 
model being considered to compute the moments of momenta. Lines beyond 
2000 must be changed by the user for each space group and model . Use of 
this section eliminates the need to enter atomic parameters for atoms 
related by symmetry to the atoms of the first molecule entered. The flow 
diagram for this program is given in figure 7. 
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Program I: COMPUTABO SUPERFICIES 


This program computes the direction cosines of the normals to 
selected crystallographic planes and calculates the areas of those 
planes . 
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100 REM PR0GRAM=C0MPUTAB0 SUPERFICIES by F'.G. Barber, <c) 1984 
110 REM Keysville, Virginia 23947 
120 REM 

130 DIM H(3) ,K(3> ,L(3) ,R<3,3) ,RP(3) 

140 REM **-*******data entry 

150 PI=3. 14159/180 

160 INPUT "Do you wish the results printed (Y/N) : " ; Y$ 

170 INPUT "Enter unit cell lengths, a,b,2<c: ";A,B,C 

180 INPUT "Enter unit cell angles, alpha, beta, S< gamma: " ; AA, BA, CA 

190 IF Y*="Y", THEN LPRINT "Unit Cell Parameters = "A; B; C; AA; BA; CA 
200 AA=AA*P I : BA=BA*F' I : CA=CA*F' I 

210 INPUT "Do you wish to enter coordinates of the points (Y/N): ";Z$ 

220 IF Z$<>"Y" THEN 270 

230 INPUT "Enter coords of first point, hi, k 1,11: " ; H ( 1 ) , K < 1 ) , L ( 1 ) 

240 INPUT "Enter coords of second point, h2, k2, 1 2: " ; H (2) , K (2) , L (2) 

250 INPUT "Enter coords of third point, h3,k3,i3: " ; H (3) , K (3) , L (3) 

260 GOTO 410 

270 INPUT "Enter Miller indices Ehkl 1 of plane: ";HX,KX,LX 

280 IF Y$="Y" THEN LPRINT "********************************** ": LPRINT " fhk 

i : = "hx;kx;lx 

290 FOR 1=1 TO 3 

300 H < I ) =0 ! : K ( I ) =0 ! : L < I ) =0 ! : NEXT I 
310 IF LX< >0 THEN 370 
320 IF HXOO THEN 350 

330 IF KX=0 THEN PRINT "ERROR)) Chkl D all zero!!":G0T0 270 ELSE K ( 1 ) =1 /KX .* H (2) =1 
: K (2) =1/KX: K (3) =1/KX: L (3) =1 
340 GOTO 410 

350 IF KXOO THEN H < 1 ) = 1 /HX : K (2) =1 /KX : H (3) =1 /HX : L (3) =1 ELSE H < 1 ) =1 /HX : H (2) =1 /HX : 

k<2)=i:h(3)=i/hx:l <3) =1 

360 GOTO 410 

370 IF HX=0 GOTO 400 

380 IF KX=0 THEN H ( 1 ) = 1 /HX : H < 2) = 1 /HX : K (2) =1 : L (3) =1 /LX ELSE H < 1 ) =1 /HX Z K <2> =1 /KX : L 
(3) =1/LX 
390 GOTO 410 

400 IF KXOO THEN H < 1 ) =1 : K < 1 ) =1 /KX: K (2) =1 /KX : L (3) =1 /LX ELSE L < 1 J =1 /LX: K (2) =1 : L (2 

)=l/LX:H(3)=l:L(3)=l/LX 

410 FOR 1=1 TO 3 

420 REM LPRINT H < I ) , K ( I ) , L < I ) 

430 H ( I ) =H ( I ) *A : K ( I ) =K ( I ) *B*S I N < CA ) : L ( I > =L < I ) *C 
440 AT=L ( I ) *C0S (BA) 

450 BT=L ( I ) *C0S ( AA) 

460 CPS=AT*AT+BT*BT+2*AT*BT*C0S (CA) 

470 L =SQR < L ( I ) *L ( I ) —CPS ) 

480 IF L ( I > < 0 THEN L < I ) =— L ELSE L ( I ) =L 
490 REM LPRINT H(I),K(I),L(I) 

500 NEXT I 

510 REM *********cal cul ate vectors 

520 FOR 1=2 TO 3 

530 R < 1 , 1 > =H < I ) -H < 1 ) : R < 1 , 2) =K < I ) -K < 1 ) : R ( 1 , 3) =L ( I ) -L < 1 > 

540 NEXT I 

550 R < 1 , 1 ) =R (3, 1) :R(1,2>=R<3,2) :.R<1,3) =R (3,3) 

560 REM *********A x B=area 

570 RP ( 1 ) =R <2, 2) *R (1,3) -R < 1 , 2) *R <2, 3) 

580 RP (2) =R ( 1 , 1 ) *R (2,3) -R (2, 1)*R(1,3> 

590 RP (3) =R (2,1) *R (1,2) — R (1,1) *R (2,2) 

600 IF HX=0 OR IF KX=0 OR IF LX=0 THEN 620 

610 FOR 1=1 TO 3:RP(I)=RP(I)/2:NEXT I 

620 AR=SQR(RP(1)*RP(1)+RP(2)*RP(2)+RP(3)*RP(3) ) 

630 REM *********print results 

640 IF Z$< >"Y" THEN 680 

650 INPUT "Is superficies triangular (Y/N): ";X$ 
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660 IF X$=" Y" , THEN AR=.5*AR 
670 GOTO 720 
680 AR=. 5*AR 

690 IF HX=0 THEN AR=2*AR : GOTO 720 
700 IF KX=0 THEN AR=2*AR : GOTO 720 
710 IF LX=0 THEN AR=2*AR 
720 PRINT 

730 PRINT "Area of face bounded by points: " 

740 FOR 1=1 TO 3 

750 PRINT H < I ) , K ( I ) , L < I ) 

760 NEXT I 

770 PRINT " is "AR" square Angstroms" 

780 PRINT 

790 PRINT "The vector perpendi cul ar to the face is " 

800 FOR 1 = 1 TO 3: RP < I ) =RP ( I ) /AR: NEXT I 

810 PRINT " N= "RP < 1 ) " i + "RP(2)"j + "RP<3)"k" 

820 PRINT 

830 IF Y$0"Y", GOTO 900 

340 LPRINT :LPRINT "Area of face bounded by points: " 

850 FOR 1=1 TO 3: LPRINT H < I ) , K < I ) , L < I > : NEXT I 

860 LPRINT ” is "AR" square Angstroms." 

870 LPRINT : LPRINT "The vector perpendi cul ar to the face is " 
880 LPRINT " N= "RP(l)"i + "RF’(2)"j + "RP(3)"k" 

890 LPRINT 

900 INPUT "Another superficiei calculation (Y/N) : ";X$ 

910 IF X$="Y", THEN 210 
920 END 
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Program II: 


MOMENTA MOMENTORUM 


This program computes the moments of momentum for a collection 
of atoms about normals to selected crystallographic planes. 
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100 REM PROGRAM = MOMENTA MOMENTORUM by P. G. Barber 
110 REM (c) Keysville, VA 23947, March 19S4 

120 REM ************************ Set-up and input data, limited to <=200 atoms** 


130 ND=200 

140 DIM H(ND) , K ( ND ) ,L(ND) , D <ND) , M (ND) ,MU<ND) ,R1 <ND) ,R2(ND> ,R3<ND) 

150 DIM RR(ND) ,RN(ND) , IN(ND) , AW(ND) ,AR(ND> , S*<ND> ,A*<ND) 

160 PI=180/3. 14159 

170 PRINT "PROGRAM: Momenta Inertiae by P.G. Barber, <0)1984” 

180 INPUT "Enter Y if you wish to print the results: " jYS 
190 INPUT "Enter D if you wish disc input of structure: "JYIS 
200 IF Y1$="D" , THEN GOSUB 1150 
210 INPUT "Enter title";TS 

220 INPUT "Enter unit cell lengths, a, b , S<c : " 5 A, B, C 

230 INPUT "Enter unit cell angles, alpha, beta, ?<gamma: " j AA, BB, CC 

240 INPUT "Enter number of atoms: ";N 

250 INPUT "Enter number of types of atoms:"; NT 

260 FOR 1=1 TO NT 

270 PRINT "For atom type # "I: INPUT " Enter atomic symbol , atomi c weight, and ato 
mic radius" 5 S$ ( I ) , AW < I ) , AR ( I ) : NEXT I 
280 CLS: PRINT: PRINT 

290 PRINT "################################################" 

300 PRINT USING " \ \";T* 

310 PRINT "################################################" 

320 PRINT " A= "A" B= "B" 0= "C 

330 PRINT "alpha= "AA" beta= "BB" gamma= "CC 

340 PRINT "Number of atoms= "N 

350 FOR 1=1 TO NT: PRINT "Atom # "I" is "AS$<I>" with wt ="AW<I>" and radius ="AR 
(I): NEXT I 

360 INPUT "Are these values correct (Y/N)?";X$ 

370 IF XSO"Y", GOTO 210 
380 IF Y$="Y" , THEN GOSUB 940 
390 AA=AA/PI : BA=BB/PI : CA=CC/PI 
400 IF Y 1 $="D" , THEN GOSUB 1400 
410 FOR 1=1 TO N 

420 PRINT "Enter ChkI 3 for atom #"Ij: INPUT HT , KT , LT 
430 H < I ) =HT*A: K < I ) =KT*B: L ( I ) =LT*C 

440 PRINT "Enter type number for atom # "I;: INPUT J 

450 PRINT "Enter multiplier of atom #" I;: INPUT MU < I ) 

460 M(I)=AW<J> :D(I)=AR(J) : A* < I > =S$ < J > 

470 PRINT : PRINT 

480 PRINT " *** atom #"I" is "A$<I)" ***" 

490 PRINT " ChkI 3= "HT, KT, LT 
500 PRINT "mass= "M(I>" g/mole" 

510 PRINT "multipl ier= "MU(I> 

520 PRINT "atomic radius= "D<I>" Angstroms" 

530 PRINT 

540 INPUT "Are these values correct (Y/N)?" ; X$ 

550 IF X$<>"Y" GOTO 420 

560 IF Y$="Y", THEN GOSUB 1050 

570 NEXT I 

580 INPUT "Enter Y if you wish to store data as file: " ; Y2$ 

590 IF Y2*="Y", THEN GOSUB 1400 
600 INPUT "Enter the ChkI 3- vector :"; Q1 , Q2, 03 
610 INPUT "Enter ChkI 3 for or i gi n : " 5 HT, KT, LT 
620 H <0> =HT*A: K <0) =KT*B: L <0) =LT*C 

630 INPUT "Enter ChkI 3 for termi num: " ; HT2, KT2, LT2 

640 H2=HT2*A: K2=KT2*B: L2=LT2*C 

650 PRINT "For vector ChkI 3 = "501,02,03 

660 PRINT "Axis is from ChkI 3= "HT, KT, LT 

670 PRINT " to ChkI 3= "HT2, KT2, LT2 

680 INPUT "Are these coordinates correct <Y/N)?";X$ 
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690 IF X$< >"Y" , GOTO 610 

700 IF Y$="Y" , THEN GOSUB 1010 

710 REM ************************ Calculation *********************************** 


720 I T=0 ! : I P 1 =0 : I P2=0 : I P3=0 

730 R 1 (O) =H2-H (O) : R2 C 0 ) =K2-K < 0 ) : R3 < 0 ) =L2-L < 0 ) 

/4o DR— SQR (R1 (u) *R1 (0) +R2 (0) *R2 (0) +R3 (O) *R3 CO) +2*R1 (0) *R2 CO) *COS CCA) +'?*Ri CO) *r= { 
0) *C0S CBA) +2*R2 CO) *R3 (0) *C0S CAA) ) ~ R1 CO> * Ro ( 

750 R1 CO) =R1 CO) /DR: R2 CO) =R2 CO) /DR: R3 CO) =R3 (O) /DR 
760 FOR 1=1 TO N 

770 R1 C I ) =H < I ) -H CO) : R2 < I ) =K C I > -K CO) : R3 C I ) =L C I ) -L CO) 

Sr, 1 SSScosIba,' 1 ’ * R3 1 1 > * R3 < 1 ’ +2 * R1 m .R3< , 

S » AA ) 


800 IN C I ) =MU C I ) *M C I ) * C . 4*D C I ) *D C I ) +RR C I ) -RN C I ) *RN CD) 

BIO IT=IT+IN(I) 

320 P 1 =MU C I ) *M C I ) *R 1 C I ) *R2 Cl): P2=MU C I ) *M C I ) *R2 ( I ) *R3 C I ) : P3=MU C I ) *M C I ) *R1 ( I ) *R3 ( I 

830 I P 1 = I P 1 +P 1 : I P2= I P2+P2 : I P3= I P3+P3 
B40 REM LPRINT I , RR < I ) , RN < I > , PI , P2, P3 
850 REM LPRINT R1 C I ) , R2 C I ) , R3 C I ) , DR 
860 NEXT I 

870 PRINT T$," 1= " ; IT 

880 PRINT "P C>:y ) = "jIPlj "PCyz) = ";IP2; "PC>:z) = " ;IP3 

890 IF Y$=" Y" , THEN GOSUB 1110 

900 INPUT "Do you wish another calculation of I CY/N):";Z1$ 

910 IF Z 1*="Y" THEN GOTO 600 
920 GOTO 1460 

930 REM ****************************** p rint 0ut initial Input Data ************ 


940 LPRINT "###############################################■■ 

950 LPRINT USING " \ \ " ; T* 

960 LPR I NT " ############################################### •• 

970 LPRINT "A= "A" B= "B" C= "C 

980 LPRINT "alpha= "AA" beta= "BB" gamma= "CC 

990 FOR 1=1 TO NT:LPRINT "Atom # "I" is "AS$CI)" with wt ="AWCI)" and radius ="A 
R Cl): NEXT I 
1000 RETURN 

1010 LPRINT : LPRINT " XXXXXXX MOMENT OF INERTIA ABOUT Chkl >="Q1 ; Q25 Q3" XXXX 

XXX" 

1020 LPRINT USING "Axis from Chkl 3=+#. #### +#.#### +#.#### to Ehkl 3=+# . #### +#.# 
### +#.####"; HT , KT , LT, HT2, KT2, LT2 
1030 LPRINT " is" 

1040 RETURN 

1050 LPRINT " **** atom #"I" ****" 

1060 LPRINT "Chkl 1= "HT, KT, LT 
1070 LPRINT "mass= "MCI)" g/mole" 

1080 LPRINT "multiplier= "MUCD 

1090 LPRINT "atomic radius= "DCI)" Angstroms" 

1100 RETURN 

1110 LPRINT ;;; "I = "IT:LPRINT "P C>:y ) =" : IP1 ; "P Cyz ) =" ; IP2* "P C>:z ) =" ; IP3 

1120 REM LPRINT "I of atom #"I" = "INC I) 

1130 RETURN 

1140 REM *********************** Disc Input of Data **************************** 


1150 INPUT "Enter data file name: " ; F$ 

1160 OPEN F$ FOR INPUT AS #1 

1170 INPUT #1, T*: INPUT #1 , A, B, C, AA, BA, CA, N, NT 

1180 FOR 1=1 TO N: INPUT #1 , H C I ) , K C I ) , L C I ) , M C I ) , MU C I ) , D ( I ) , A$ C I ) : NEXT I 
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1190 REM *******************To Subroutine to expand with sym- related atoms**** 

1200 REM GOSUB 2000 

1210 IF Y$< >" Y" , THEN GOTO 1370 

1220 REM *********************** Disc Input Print Out ************************** 
1230 LPRIMT : LPRINT 

1240 LPRINT USING " \ \ n ;T^ 

1250 LPRINT : LPRINT 

1260 LPRINT "Disc entry of data -from data -file #: "F$ 

1270 LPRINT "a = "A"b = M B"c = "C 

1280 LFRINT "alpha = " (AA*PI) " beta = "<BA*PI>” gamma = "(CA*PI> 

1290 LPRINT N" atoms entered in structure. " 

1300 LPRINT "*** atom # *** Ehkl 1 mass mu D(A>" 

1310 V*=" ### \ \ +##.#### +##.#### +##.#### ###.### #.## #.##" 

1320 INPUT "Enter multiplier factor:"?P 
1330 FOR 1=1 TO N: MU < I ) =MU ( I ) *P2 NEXT I 
1340 FOR 1=1 TO N 

1350 LPRINT USING VS 5 I A* ( I ) , H < I ) /A, K < I ) /B, L < I ) /C, M < I ) , MU < I ) , D ( I ) 

1360 NEXT I 
1370 RETURN 600 

1380 FOR 1=1 TO N:MU<I)=MU<I)*P:NEXT I 
1390 RETURN 610 

1400 INPUT "Enter data disc 'File name: ";F1$ 

1410 OPEN FI* FOR OUTPUT AS #2 

1420 PRINT #2, T*: PRINT #2, A, B, C, AA ? BA, CA„ N, NT 

1430 FOR 1=1 TO N: PRINT #2 P H < I > , K < I ) , L ( I ) , M < I ) , MU < I ) , D ( I ) , A* < I > : NEXT I 
1440 CL0SE#2 
1450 RETURN 
1460 END 

2000 REM *******************Subrouti ne expands input data set******************* 

2010 REM *******************Change as needed for space group or model ********** 
2020 H (N+l >=H < 1 ) +A/2:K (N+l ) =K (1 ) : L (N+l ) =L ( 1 ) +C/2 : H ( N+2 ) =H < 2 ) +A/2: K <N+2) =K <2) : L <N 
+2) =L (2) “C/2 
2030 N=2*N 

2040 H <N+1 ) =H < 1 ) :K(N+1)=K(1) : L <N+1 > =L < 1 ) +C: H <N+2) =H (2) +A/2: K <N+2) =K <2) :L<N+2)=L< 
2) +C/2 
2050 N=N+2 

2060 H(N+1)=H<1)+A:K(N+1)=K(1) :L(N+1)=L(1)+C:H<N+2>=H<2)+A:K<N+2)=K<2) : L (N+2) =L ( 
2) : H (N+3) =H < 1 ) +A: K <N+3> «K ( 1 ) : L <N+3> =L < 1 ) 

2070 N=N+3 

2080 FOR 1=10 TO 18 

2090 H < I ) =H ( I -9 > : K ( I ) =K ( I -9 ) +B /2 : L ( I ) =L ( I -9 > 

2100 M < I ) =M ( 1—8) :MU<I)=MU(I-8) :D(I)=D(I-8) : A$(I)=A$<I-8> 

2110 NEXT I 

2120 N=N+9 

2130 FOR 1=19 TO 27 

2140 H < I ) =H <1 — 18) :K<I) =K (1 — 18) +B: L ( I ) =L (1 — 18) :M(I) =M (1 — 18) : MU < I ) =MU (1-18) : D ( I > =D 
(1—18) :A$(I)=A$(I-18) 

2150 NEXT I 
2160 N=N+9 

2170 FOR 1=3 TO 27 STEP 2 

2180 M ( I ) =M ( 1 ) :MU(I)=MU<1) : D ( I ) =D ( 1 ) : A$(I)=A*(1) 

2190 M ( I+i ) =M (2) :MU(I+1)=MU(2> :D(I+1)=D(2) : A$ < 1+1 ) =A* (2) 

2200 NEXT I 
2210 GOTO 2330 

2220 REM ************************** Subroutine to select ions***************** 
2230 H(1>=H<3> :K<1)=K<3> : L ( 1 ) =L (3) : M ( 1 ) =M (3) :MU(1)=MU(3) :D(1)=D(3) :A$(1)=A$(3> 
2240 H(2)=H<4> :K(2)=K<4> :L(2)=L<4) :M<2)=M<4) : MU (2) =MU (4) :D(2)=D<3) :A$(2)=A$<4> 
2250 H(3)=H(10):K(3)=K(10):L(3)=L(10):M(3)=mi0) : MU (3) =MU ( 10) : D <3) =D (10) : A$ (3) =A 
$(10) 
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2260 H < 4 ) =H ( 1 1 > : K ( 4 ) =K < 1 1 ) : L < 4 ) =L < 1 1 ) : M ( 4 ) =M ( 1 1 ) : MU ( 4 ) =MU < 1 1 ) : D < 4 ) =D < 1 1 ) : A$ < 4 ) =A 
$(11) 

2270 H <5) =H < 16) : K (5) =K < 16) :L(5)=L<16> : M (5) =M ( 16) : MU (5) =MU < 16) : D (5) =D < 16) : A$ (5) =A 
$ (16) 

2280 H<6)=H(17):K(6) =K (17) :L(6) =L (17) :M(6) =M (17) .* MU <6 ) =MU ( 1 7) : D (6) =D < 17) : A$ (6) =A 
$(17) 

2290 H(7)=H(21) :K(7)=K(21) :L(7)=L<21) :M(7)=M<21) : MU (7) =MU (21 ) : D (7) =D (21 ) : A$ (7) =A 
$(21) 

2300 H (8) =H (24) : K (8) =K (24) : L (8) =L (24) : M (8) =M (24) : MU (8) =MU (24) : D (8) =D (24) i A$ (8) =A 
$ (24) 

2310 H ( 9 > =H (12) : K ( 9 ) =K (12) : L ( 9 ) =L (12) :M(9) =M (12) *.MU (9) =MU( 12) : D (9) =D ( 12) : A$ (9) =A 
$(12) 

2320 N=9 
2330 RETURN 


31 



REFERENCES 


1. J. D. H. Donnay; and D. Harker: "A New Law of Crystal Morphology 

Extending the Law of Bravais," Am. Mineralologist 22 (1937) 446. 

2. J. W. Gibbs: Collected Works (Longmans, Green and Company, New York) 

(1928). Footnote page 325. 

3. G. Wul ff : Z. Krist. 34 (1901) 449. 

4. G. H. Gilmer; and P. Bennema: "Computer Simulation of Crystal Surface 

Structure and Growth Kinetics," J. Crystal Growth, 13/14 (1972) 148. 

5. G. H. Gilmer: "Computer Simulation of Crystal Growth," J. Crystal 

Growth 42 (1977) 3. 

6. P. P. von Veimarn: "The Precipitation Laws," Chem. Rev. 2 (1925) 217. 

7. J. M. Bermond; and J. A. Venables: "A Model for the Growth Shape of 
Crystallites on Surfaces," J. Crystal Growth 64 (1983) 239. 

8. I. N. Stranski; and R. Kaischev: "Equilibrium Forms of Homeopolar 

Crystals," Z. Krist. 78 (1931) 373. 

9. I. N. Stranski; and R. Kaischev: "Crystal Growth and Formation of 

Nuclei," Physik. Z. 36 (1935) 393. 

10. W. K. Burton; and N. Cabrera: "Crystal Growth and Surface Structure," 

Disc. Faraday Soc. 5 (1949) 33. 

11. W. K. Burton; N. Cabrera; and F. C. Frank: "Growth of Crystals and 

the Equilibrium Structure of Their Surfaces," Phil. Trans. Roy. Soc. 
(London) 243A (1951) 299. 

12. H. E. Buckley: Crystal Growth , John Wiley & Sons, N.Y. (1951). 

13. P. Hartman: "Equilibrium Forms of Some Ionic Crystals," Neues Jahrb. 

Mineral. Monatsh. (1959) 73. 

14. G. A. Wolff; and J. G. Gualtieri: "PBC (Periodic Bond Chain) Vector, 

Critical Bond Energy Ratio, and Crystal Equilibrium Form." Am. 
Mineralogist 47 (1962) 562. 

15. W. Kleber: "Structure Relations and Principal Bond Vectors," 

Tonind.-Ztg. Keram. Rundsch 86 (1962) 251. 

16. H. Honigmann: "Equilibrium Forms and Growth Forms of Crystals," 

Bull. Soc. Franc. Mineral Crist. 86 (1963) 408. 

17. H. Hilton: Mathematical Crystallography , Dover Publishers, N.Y., 

(1963), Chapter XI. 


32 



18. A. I. Kitaigorodsky: Molecular Crystals and Molecules, Academic 

Press, N.Y. (1973). 

19. P. Hartman; and W. G. Perdok: "On the Relations Between Structure and 

Morphology of Crystals I." Acta. Cryst. 8 (1955) 49. 

20. P. Hartman; and P. Bennema: "The Attachment Energy as a Habit 

Controlling Factor," J. Crystal Growth 49 (1980) 145. 

21. D. Aquilano; M. Franchini -Angela; M. Rubbo; G. Mantovani ; and G. 

Vaccani : "Growth Morphology of Polar Crystals - A Comparison Between 

Theory and Experiment in Sucrose," J. Crystal Growth 61 (1983) 369. 

22. Michael Saska; and Allan S. Myerson: "The Theoretical Shape of 

Sucrose Crystals from Energy Calculations," J. Crystal Growth 61 
(1983) 546. 

23. R. Hamar; and S. Hamar-Thibault: "Observations sur la Croissance du 

Compose CrB," J. Crystal Growth 63 (1983) 365. 

24. R. J. Davey; J. W. Mullin; and M. J. L. Whiting: "Habit Modification 

of Succinic Acid Crystals Grown from Different Solvents," J. Crystal 
Growth 58 (1982) 304. 

25. H. J. Human; J. P. van der Eerden; L. A. M. J. Jetten; and J. G. M. 

Odekerken: "On the Roughening Transition of Biphenyl: Transition of 

Faceted to Non-Faceted Growth of Biphenyl for Growth from Different 
Organic Solvents and the Melt," J. Crystal Growth 51 (1981) 589. 

26. V. V. Podlinsky; and V. G. Drykin: "Surface Structure Change of a 

Crystal in Solution," J. Crystal Growth 62 (1983) 532. 

27. C. S. Strom; and W. M. M. Heijnen: "Computerized PBC Analysis of 

Weddellite," J. Crystal Growth 51 (1981) 534. 

28. J. J. Burton: "Configuration, Energy, and Heat Capacity of Small 

Spherical Clusters of Atoms," J. Chem. Phys. 52 (1970) 345. 

29. J. J. Burton: "Densely Packed Small Clusters of Atoms," Nature 229 

(1971) 335. 

30. M. Hoare; and P. Pal: "Physical Cluster Mechanics. Statics and 

Energy Surfaces for Monatomic Systems." Ad van. Phys. 20 (1971) 161. 

31. M. R. Hoare; and P. Pal: Statics and Stability of Small Cluster 

Nuclei ," Nature (London), Phys. Sci . , 230 (1971) 5. 

and 

M. R. Hoare; and P. Pal: "Geometry and Stability of Spherical fee. 

Microcrystallites, 11 Nature (London), Phys. Sci., 236 (1972) 35. 


33 



32. M. Hoare; and P. Pal: "Statistics and Stability of Small Assemblies 

of Atoms," J. Crystal Growth 17 (1972) 77. 

33. D. Bonchev; 0. Mekenyan; and H. Fritsche: "An Approach to the 

Topological Modeling of Crystal Growth," J. Crystal Growth 49 (1980) 
90. 


34. C. Kittel : Intro, to Solid State Physics , Third Edition, John Wiley 

and Sons, Inc., N.Y. (1966). 

35. A. Julg; and B. Deprick: "The Influence of Glycine, According to pH, 
Upon the Growth Habit of NaCl and KI, in Aqueous Solution," J. Crystal 
Growth 62 (1983) 587. 

36. J. D. Birchall ; and R. J. Davey: "The Crystallization of Sodium 

Chloride from Aqueous Solution in the Presence of Polysaccharides," 

J. Crystal Growth 54 (1981) 323. 

37. J. S. Broadley; D. W. J. Cruickshank; J. D. Morrison; J. M. Robertson 

(FRS); and H. M. M. Shearer: "Three-Dimensional Refinement of the 

Structure of -Succininc Acid," Proc. Roy. Soc. (London) 251A (1959) 
441. 

38. G. M. Brown; and H. A. Levy: "Further Refinement of the Structure of 

Sucrose Based on Neutron-Diffraction Data," Acta. Cryst. B29 (1973) 
790. 


39. A. F. Wells: "Crystal Habit and Internal Structure. I, II, and III." 

Phil. Mag. 37 (series #7) (1946) 184, 217, 605. 

40. Grant R. Fowles: Analytical Mechanics. Holt, Rinehart and Winston, 

Inc. (1962). 


34 



Table 1 

Summary of Calculations for Sodium Chloride 
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Table 2 

Summary of Calculations for -Succinic Acid 
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Table 3 

Summary of Calculations for Sucrose 
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Table 4 

Summary of Calculations on Chromium Boride 
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Figure 6 1 flowchart for progror* Cornputobo Superficies/ poge 1 
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Figure 6s flowchart for progrom Cornputobo Superficies# poge 2 
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Figure 7 1 flowchart for progron Momento Monentorum# page 1 
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